ON COMPUTATIONAL COMPLEXITY OF SIEGEL JULIA SETS 
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Abstract. It has been previously shown by two of the authors that some polynomial 
Julia sets are algorithmically impossible to draw with arbitrary magnification. On the 
other hand, for a large class of examples the problem of drawing a picture has polynomial 
complexity. In this paper we demonstrate the existence of computable quadratic Julia 
sets whose computational complexity is arbitrarily high. 



1. Foreword 

Let us informally say that a compact set in the plane is computable if one can program 
a computer to draw a picture of this set on the screen, with an arbitrary desired magnifica- 
tion. It was recently shown by the second and third authors, that some Julia sets are not 
computable [BYJ. This in itself is quite surprising to dynamicists - Julia sets are among 
the "most drawn" objects in contemporary mathematics, and numerous algorithms exist 
to produce their pictures. In the cases when one has not been able to produce informative 
pictures (the dynamically pathological cases, like maps with a Cremer or a highly Liou- 
ville Siegel point) the feeling had been that this was due to the immense computational 
resources required by the known algorithms. 

The next surprise came with the discovery by the authors of this paper in BBY that 
all Cremer quadratics (or more generally, rational maps without rotation domains) have 
computable Julia sets. The non-computable examples constructed in [BYj were Siegel 
quadratic polynomials, and one would expect the Cremer case to be at least as bad if not 
worse computationally. 

The natural question to ask is then whether in those cases in which we know the Julia 
set is computable, but no good pictures exist, the computational complexity of such a set 
is indeed high. Here at least, our original intuition seems to be correct: it is shown in the 
present paper that there exist computable Siegel quadratic Julia sets with arbitrarily high 
computational complexity An irritating possibility still remains that some Cremer Julia 
sets are computationally easy (and we just do not go about trying to draw them in the right 
way). This, however, seems unlikely. We note that the examples constructed in this paper 
are the first known cases of Julia sets which are not poly-time computable. The second 
author |Brvlj and independently Rettinger |Retj have previously shown that hyperbolic 
Julia sets are poly-time computable. More recently the second author has shown |Brv2j 
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that some Julia sets with parabolics are poly-time computable as well. The last result was 
yet another surprise, as the time complexity of all previously known algorithms for these 
Julia sets was exponential. 

The structure of the paper is as follows. In ^2.21 of the Introduction, having stated the 
principal definitions, we formulate the main result of the paper. In $ 12.41 we give a sketch 
of the argument. In §3] we prove several technical lemmas. The final ^contains the proof 
of the Main Theorem. 

Acknowledgement. We would like to thank Giovanni Gallavotti for very helpful sugges- 
tions on the exposition. 

2. Introduction 

2.1. Computability of real sets. The reader is directed to |BYj for a more detailed 
discussion of the notion of computability of subsets of M n as applied, in particular, to Julia 
sets. We recall the principal definitions here. The exposition below uses the concept of a 
Turing Machine. This is a standard model for a computer program employed by computer 
scientists. Readers unfamiliar with this concept should think instead of an algorithm 
written in their favorite programming language. These concepts are known to be equivalent. 

Denote by D the set of the dyadic rationals, that is, rationals of the form We say 
that : N — > D is an oracle for a real number x, if \x — 4>{n)\ < 2~ n for all n G N. In 
other words, provides a good dyadic approximation for x. We say that a Turing Machine 
(further abbreviated as TM) is an oracle machine, if at every step of the computation 
M is allowed to query the value <fi{n) for any n. This definition allows us to define the 
computability of real functions on compact sets. 

Definition 2.1. We say that a function / : [a, b] — > [c,d] is computable, if there exists an 
oracle TM M^(m) such that if <fi is an oracle for x G [a, b], then on input m, outputs 
at/6D such that \y — f(x)\ < 2~ m . 

To understand this definition better, the reader without a Computer Science background 
should think of a computer program with an instruction 

READ real number x WITH PRECISION n(m). 

On the execution of this command, a dyadic rational d is input from the keyboard. This 
number must not differ from x by more than 2~ nt - m ' (but otherwise can be arbitrary). The 
algorithm then outputs f(x) to precision 2~ n . 

It is worthwhile to note why the oracle mechanism is introduced. There are only countably 
many possible algorithms, and consequently only countably many computable real numbers 
which such algorithms can encode. Therefore, one wants to separate the hardness of 
encoding the real number x from the hardness of computing the value of the function f(x), 
having the access to the value of x. 

Let K C M fe be a compact set. We say that a TM M computes the set K if it approximates 
K in the Hausdorff metric. Recall that the Hausdorff metric is a metric on compact subsets 
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of E n defined by 
(2.1) 



d H (X,Y) 



inf{e > | X C U e (Y) and Y C U e (X)}, 



where U e (S) is defined as the union of the set of e-balls with centers in S. 
We introduce a class C of sets which is dense in metric du among the compact sets and 
which has a natural correspondence to binary strings. Namely C is the set of finite unions 
of dyadic balls: 



Members of C can be encoded as binary strings in a natural way. 

We now define the notion of computability of subsets of M. n (see |Weij . and also |RWj ). 

Definition 2.2. We say that a compact set K C M fc is computable, if there exists a TM 
M(d,n), where d G D, n G N which outputs a value 1 if dist(<i, K) < 2 _n , the value if 
dist(c/, K) > 2 • 2~ n , and in the "in-between" case it halts and outputs either or 1. 

In other words, it computes, in the classical sense, a function from the family Tk of 
functions of the form 



Theorem 2.1. For a compact K C M fc the following are equivalent: 

(1) K is computable as per definition \2.2i 

(2) there exists a TM Mim), such that on input m, M(m) outputs an encoding ofC m G C 
such that dniKjCm) < 2~ m (global computability), 

(3) the distance function dx{x) = inf{|x — y\ \ y G K} is computable as per definition 



Note that in the case k = 2 computability means that K can be drawn on a computer 
screen with arbitrarily good precision (if we imagine the screen as a lattice of pixels). 

In the present paper we are interested in questions concerning the computability of the 
Julia set J c = J(f c ) = J(z 2 + c). Since there are uncountably many possible parameter 
values for c, we cannot expect for each c to have a machine M such that M computes J c 
(recall that there are countably many TMs). On the other hand, it is reasonable to want 
M to compute J c with an oracle access to c. Define the function J : C — > K* (K* is the 
set of all compact subsets of C) by J(c) = J(f c )- In a complete analogy to Definition 12.11 
we can define 

Definition 2.3. We say that a function k : S — > K* for some bounded set S is computable, 
if there exits an oracle TM M^(d,n), where <fi is an oracle for x G S, which computes a 
function ()2.2|) of the family T K {x). 

Equivalently, there exists an oracle TM M^m) with <fi again representing x G S such 
that on input m, outputs aCeC such that d H (C, n(x)) < 2~ m . 




(2.2) 




ro 
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In the case of Julia sets: 

Definition 2.4. We say that J c is computable if the function J : d h- > is computable 
on the set {c}. 



We have the following (see |Brvlj ): 



Theorem 2.2. Suppose that a TM computes the function J on a set S C C. T/ien J 
zs continuous on S in Hausdorff sense. 

Proof. Let c be any point in S, and let e = 2~ k be given. Let <fi be an oracle for c such 
that \4>(ri) — c| < 2-(" +1 ) for all k. We run M^ik + 1) with this oracle <p. By the definition 
of J, it outputs a set L which is a 2~( fc+1 ) approximation of J c in the Hausdorff metric. 

The computation is performed in a finite amount of time. Hence there is an m such that 
is only queried with parameters not exceeding m. Then for any x such that \x—c\ < 2~( m+1 \ 
is a valid oracle for x up to parameter value of m. In particular, we can create an oracle 
ijj for x that agrees with on 1, 2, . . . , m. If x & S, then the execution of M^ik + 1) 
will be identical to the execution of M^ik + 1), and it will output L which has to be an 
approximation of J x . Thus we have 

d H ( J c , J x ) < d H (J c , L) + d H {J x , L) < 2~( fc+1 ) + 2-( fc+1 > = 2- k . 

This is true for any x E B(c, 2~^ m+1 ^) R S. Hence J is continuous on S. □ 

The second and third authors have demonstrated in IBYI: 



Theorem 2.3. There exists a parameter value c G C such that the Julia set of the quadratic 
polynomial f c {z) = z 2 + c is not computable. 

The quadratic polynomials in Theorem 12.31 possess Siegel disks (see £ 12.31 below for the 
definitions of Siegel and Cremer points). It was further shown by the authors of the present 
paper in |BBY| that the absence of rotation domains, that is either Siegel disks or Herman 



rings, guarantees computability of the rational Julia set. This implies, in particular, that 
all Cremer quadratic Julia sets are computable - this despite the fact that no informative 
high resolution images of such sets have ever been produced. One expects, however, that 
such "bad" but still computable examples have high algorithmic complexity, which makes 
the computational cost of producing such a picture prohibitively high. We note, that the 
second author |Brvlj and independently Rettinger |Retj have shown: 

Theorem 2.4. Hyperbolic Julia sets are computable in polynomial time. That is, if J is 
the Julia set of a hyperbolic rational mapping R, then there exists a TM M(d, n) which 
computes a function of the family in time polynomial in the bit size of the input 

(d,n). It is worth noting that the same oracle TM M^{d,n) with the oracle representing 
the parameters of the rational mapping R, can be selected for all hyperbolic Julia sets of 
the same degree. Moreover, the asymptotics of the polynomial time bound depends only on 
R but not on the input (d,n). 
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2.2. Statement of the Main Theorem. On the other end of the complexity spectrum 
we expect to find "bad" but computable Siegel Julia sets and Cremer Julia sets. Indeed, 
it the present paper we show: 

Theorem 2.5. There exist quadratic Siegel Julia sets of arbitrarily high computational 
complexity. More precisely, for any computable increasing function h : N — > N there exists 
a computable Siegel parameter value c G C such that: 

• the Julia set J c is computable by an oracle TM; 

• for any oracle TM M^(m) which computes the 2~ m - approximations to J c , there 
exists a sequence {rrii\ c *L l such that requires the time of at least h(rrii) to compute 
the approximation C„ H G C. 

From this statement for global computational complexity immediately follows the corre- 
sponding local statement: 

Corollary 2.6. There exist computable parameter values c for which the Julia set J c is 
computable, and the complexity of the problem of computing a function \2.ty) in the family 
Tj c is arbitrarily high. 

2.3. Siegel disks of quadratic maps. Let R : C — > C be a rational map of the Riemann 
sphere. For a periodic point z = R p (z ) of period p its multiplier is the quantity A = 
A(^o) = DR p (zq). We may speak of the multiplier of a periodic cycle, as it is the same 
for all points in the cycle by the Chain Rule. In the case when |A| 7^ 1, the dynamics 
in a sufficiently small neighborhood of the cycle is governed by the Intermediate Value 
Theorem: when < |A| < 1, the cycle is attracting [super- attracting if A = 0), if |A| > 1 
it is repelling. Both in the attracting and repelling cases, the dynamics can be locally 
linearized: 



where ip is a conformal mapping of a small neighborhood of zq to a disk around 0. By 
a classical result of Fatou, a rational mapping has at most finitely many non-repelling 
periodic orbits. 

In the case when A = e 2me , 9 G R, the simplest to study is the parabolic case when 
9 = n/m G Q, so A is a root of unity. In this case RP is not locally linearizable; it is not 
hard to see that z G J(R). In the complementary situation, two non- vacuous possibilities 
are considered: Cremer case, when RP is not linearizable, and Siegel case, when it is. 
In the latter case, the linearizing map ip from ()2.3|) conjugates the dynamics of RP on a 
neighborhood U(zq) to the irrational rotation by angle 9 (the rotation angle) on a disk 
around the origin. The maximal such neighborhood of zq is called a Siegel disk. 

Let us discuss in more detail the occurrence of Siegel disks in the quadratic family. For 
a number 9 G [0, 1) denote [r ,ri, . . . ,r n , . . .}, r { G NU {00} its possibly finite continued 



(2.3) 
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fraction expansion: 

(2.4) [r ,n, 



r H 

ri + 



1 

■ + 

r n + ■ ■ 



Such an expansion is defined uniquely if and only if 9 ^ Q. In this case, the rational 
convergents p n /<ln = [tq, ■ ■ ■ , r n _i] are the closest rational approximants of 9 among the 
numbers with denominators not exceeding q n . In fact, setting A = e 27r , we have 

\\ h - 1| > \\ qn - 1| for all < h< q n+1 , h ^ q n . 

The difference \X qn — 1| lies between 2/g n+1 and 27i~/g n+1 , therefore the rate of growth of 
the denominators q n describes how well 9 may be approximated with rationals. 
We recall a theorem due to Brjuno (1972): 

Theorem 2.7 ( jBruj ). Let R be an analytic map with a periodic point zq G C. Suppose 
that the multiplier of z is A = e 2m0 , and 

(2.5) ^)=^^n ±1 ) <oo 

Then z is a Siegel point. 

Note that a quadratic polynomial with a fixed Siegel disk with rotation angle 9 after an 
affine change of coordinates can be written as 

(2.6) P e (z) =z 2 + e 2nie z. 



In 1987 Yoccoz [Yoq proved the following converse to Brjuno's Theorem: 

Theorem 2.8 ( |Yocj ). Suppose that for 9 e [0, 1) the polynomial Pg has a Siegel point at 
the origin. Then B(9) < oo. 

The numbers satisfying (|2.5J) are called Brjuno numbers; the set of all Brjuno numbers will 
be denoted B. It is a full measure set which contains all Diophantine rotation numbers. 
In particular, the rotation numbers [ro, ri, ...] of bounded type, that is with suprj < oo 
are in B. The sum of the series ()2.5|) is called the Brjuno function. For us a different 
characterization of B will be more useful. Inductively define 9\ = 9 and 9 n+ i = {l/9 n }. In 
this way, 

9-n ['"n— lj '"nj 5™n+l) ■ ■ ■]■ 

We define the Yoccoz 's Brjuno function as 



oo 1 

${9) = Y j 9 1 9 2 ---9 n _ l \og-. 
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One can verify that 

B{9) <oo^ $(0) < oo. 
The value of the function $ is related to the size of the Siegel disk in the following way. 

Definition 2.5. Let (U, u) be a simply-connected subdomain of C with a marked interior 
point. Consider the unique conformal isomorphism <f> : D i— > U with 0(0) = u, and 
0'(O) > 0. The conformal radius of (U, u) is the value of the derivative r(U, u) = 0'(O). 

Let P{9) be a quadratic polynomial with a Siegel disk 9 0. The conformal radius 
of the Siegel disk Ag is r{6) = r(Ag,0). For all other 9 G [0, oo) we set r{6) = 0, and 
A e = {0}. 

By the Koebe 1/4 Theorem of classical complex analysis (see e.g. |Ahlj ) . the radius of the 
largest Euclidean disk around u which can be inscribed in U is at least r(U,u)/A. 
We note that one has the following direct consequence of the Caratheodory Kernel Theorem 
(see e.g. [Pom]): 

Proposition 2.9. The conformal radius of a quadratic Siegel disk varies continuously with 
respect to the Hausdorff distance on Julia sets. 

Yoccoz |Yocj has shown that the sum 

$(0) + logr(0) 

is bounded below independently of 6 B. Recently, Buff and Cheritat have greatly 
improved this result by showing that: 

Theorem 2.10 ( BC ). The function 9 \— > $(#) + logr(#) extends to R as a 1-periodic 
continuous function. 

In |BBYj we obtain the following result on computability of quadratic Siegel disks: 

Theorem 2.11. The following statements are equivalent: 

(I) the Julia set J{Pe) is computable; 
(II) the conformal radius r(9) is computable; 
(III) the inner radius m.i z ^Q^ g \ z\ is computable. 

We note that when 9 is not a Brjuno number, the quantities in (II) and (III) are each equal 
to zero, and the claim is simply that J(Pe) is computable in this case. 

We will make use of the following Lemma which bounds the variation of the conformal 
radius under a perturbation of the domain. It is a direct consequence of the Koebe Theorem 
(see e.g. |RZj for a proof). 

Lemma 2.12. Let U be a simply- connected subdomain of C containing the point in the 
interior. Let V C U be a subdomain of U. Assume that dV C B t (dU). Then 

< r(U, 0) - r(V, 0) < Ay/r(U,0)y/e. 
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2.4. Outline of the construction. We can now describe the idea of our construction. 
This outline is rather sketchy and suffers from obvious logical deficiencies, however, it 
presents the construction in a simple to understand form. Consider the oracle Turing 
machines with representing the parameter in Pq. Since there are only countably 
many Turing machines, we may order these machines in a sequence Mf, Mf , ... We denote 
S{ the domain on which Mf computes J{Pq) properly. We thus have that for each i, the 
function J : i— > J(P$) is continuous on Si. 

Let us start with a machine M% which computes J(Pe„) for 9* = [1, 1, 1, . . .]. If any of 
the digits r« in this infinite continued fraction is changed to a sufficiently large N G N, the 
conformal radius of the Siegel disk will become small. For N — > oo the Siegel disk will 
implode and its center will become a parabolic fixed point in the Julia set (see |Do2j ). 

If we are careful, we may select %\ > 1 and N\ » 1 in such a way, that for 9\ given by 
the continued fraction where all digits are ones except r ix = Ni we have 

(2.7) r(0*)(l - 1/4) < r{9 x ) < r(0*)(l - 1/8). 

By the Koebe 1/4-Theorem, there exists l\ > such that the distance between the two 
Julia sets 

d H (J(P et ),J(P dl ))>2- e \ 

To ensure that the machine M% will not be able to produce an accurate 2~^ 1 -approximation 
of J(Pei) faster than in the time h(£i) we simply select %\ > h(£i). This guarantees that 
the TM will have to read at least h(£i) digits of the oracle <fi to distinguish the two Julia 
sets, which takes the time h(£i). 

To "fool" the machine M% 2 we then change a digit for 12 > i\ sufficiently far in the 
continued fraction of 9\ to a large A^2- In this way, we will obtain a Brjuno number 02 for 
which 

(2.8) r(0,)(l - 1/4 - 1/8) < r(0 2 ) < r(0*)(l - 1/4). 
Again, there exists £2 such that for any such Brjuno number, we have 

d H (J(Pe 1 ),J(Pe 2 ))>2- i \ 

and we choose 12 > h{£.2\ Continuing inductively, we arrive at the desired limiting Brjuno 
number 0^. 

To convince the reader that this construction is not artificial, and not due to the pecu- 
liarities of the selected computation model let us recast it somewhat informally as follows. 
It is possible by an arbitrarily small perturbation of the parameter to cause a detectable 
disturbance in the picture of J{Pe). To distinguish the picture of the new Julia set from 
the old one, in practice one needs to draw it with arbitrary precision arithmetic. That is, 
not only the input of the parameter (reading the oracle) will take a long time due to the 
number of significant digits, but also all the arithmetic manipulations with this parameter. 
Of course, the former consideration is already sufficient to prove the theorem. 
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3. Computing Noble Siegel Disks 

The primary goal of the present paper is to show that there are computationally hard 
yet computable Julia sets with Siegel disks. To establish this computability we need a 
computability result for noble Siegel disks. The term "noble" is applied in the literature 
to rotation numbers of the form [do, Oi, ■ ■ ■ , ctfc, 1, 1,1,...]. The noblest of all is the golden 
mean 7* = [1, 1, 1, . . .]. 

Lemma 3.1. There is a Turing Machine M , which given a finite sequence of numbers 
[do, Oi, . . . , dfc] computes the conformal radius r 7 for the noble number^ = [a , . . . , a^, 1, . . .] . 

The idea is to approximate the boundary of A 7 with the iterates of the critical point c 7 = 
—e 2nn /2. It is known that in this case the critical point itself is contained in the boundary. 
The renormalization theory for golden-mean Siegel disks (constructed in |McM| ) implies 
that the boundary A 7 „ is self-similar up to an exponentially small error. In particular, 
there exist constants C > and A > 1 such that 

d H ({P*ScyJ,i = 0, . . . ,U,<9A 7 J < C\~ n 

Below we derive a similar estimate for all noble Siegel disks with constructive constants 
C and A. For this, we do not need to invoke the whole power of renormalization theory. 
Rather, we will use a theorem of Douady, Ghys, Herman, and Shishikura |Dolj which 
specifically applies to quadratic noble Siegel disks. 

Noble (or more generally, bounded type) Sigel quadratic Julia sets may be constructed 
by means of quasiconformal surgery on a Blaschke product 

/ -J Z ) = e ^{ l)z 2^_ 

Jiy ! 1-32 
This map homeomorphically maps the unit circle T onto itself with a single (cubic) critical 
point at 1. The angle r(7) can be uniquely selected in such a way that the rotation number 
of the restriction p(/ 7 |t) = 7- 
For each n, the points 

{lj^.fll),..,^)} 

form the n-th dynamical partition of the unit circle. We have (cf. Theorem 3.1 of [dFdMj 
the following: 

Theorem 3.2 (Universal real a priori bound). There exists an explicit constant B > 1 
independent 0/7 and n such that the following holds. Any two adjacent intervals I and J 
of the n-th dynamical partition of f 1 are B- commensurable: 

B~ l \I\ < \J\ < B\I\. 

Let us now consider the mapping \1/ which identifies the critical orbits of / 7 and P 1 by 
We have the following (Theorem 3.10 of |YZj ) : 
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Theorem 3.3 (Douady, Ghys, Herman, Shishikura). The mapping ^ extends to a 
K -quasiconformal homeomorphism of the plane C which maps the unit disk D onto the 
Siegel disk A 7 . The constant K depends on B and do, ... , a& in a constructive fashion. 

Elementary combinatorics implies that each interval of the n-th dynamical partition 
contains at least two intervals of the (n + 2)-nd dynamical partition. This in conjunction 
with Theorem 13.21 implies that the size of an interval of the (n + 2)-nd dynamical partition 
of f 1 is at most r n where 

B 



B + l 

We now complete the proof of Lemma 13.11 Denote W n the connected component con- 
taining of the domain obtained by removing from the plane a closed disk of radius 2Kr n 
around each point of 

= {P^(Cy), 2 = 0,..., q n+2 }- 

By Theorem 13. 3| 

dist H (fi n ,<9A 7 ) < Kr n , 

and we have 

W n C A 7 and dist H (dA 1} dW n ) < e n = 2Kr n . 

Any constructive algorithm for producing the Riemann mapping of a planar region (e.g. 
that of u BBj) can be used to estimate the conformal radius r(W n , 0) with precision e n . 
Denote this estimate r n . 

Elementary estimates imply that the Julia set J(Pj) C -8(0, 2). By Schwarz Lemma this 
implies r(A 7 ,0) < 2. By Lemma f2. 121 we have 

|r(A 7 ,0)-r„| < \r(A 7 ,0)-r(W n ,0)\ + \r(W n ,0)-r n \ < A^T n + e n — ■+ 0, 
and the proof is complete. 

4. Making Small Changes to $ and r 
For a number 7 = [a 1; a 2 , . . .] G K \ Q we denote 

oii{l) = , 



O'i+l + 



&i+2 + • • • 

so that 



Hi) = 53 ai (7)^2(7) • • -0^-1(7) !og^r- 



n>l 

We will show the following two lemmas. 
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Lemma 4.1. For any initial segment I — ycLo, di, . . . , a n ), write lu = [do, Oi, • • • , d n , 1,1,1,...]. 
Then for any e > 0, there is an m > and an integer N such that if we write ft = 
[do, at, ... , a n , 1, 1, . . . , 1, N, 1, 1, . . .], where the N is located in the n + m-i/j position, then 

Lemma 4.2. Foru as above, for any e > inere is an mo > 0, which can be computed from 
(a , ai, . . . , a n ) and e, st/cn i/iai /or any m > m , and for any tail I = [a n+m , a n+m+ x, . . .] 
if we denote 

ft [^1; ® , 2i ■ ■ ■ i Q"ni 1) 1, • • • , 1) <^n+m, Q"n+m+li ■ ■ ■] , 

then 

$(/3 J ) > -e. 
We first prove lemma UTTl Denote 

$-(^) = - a (^)«i(^) • ..an+m-i(w) log 

The value of the integer m > is yet to be determined. Denote 

ft N = (a , di, . . . , a n , 1, 1, . . . , 1, AT, 1,1,...). 
We will need the following estimates, which are proven by induction 

Lemma 4.3. For any N, the following holds: 
(1) For i < n + m we have 



(2) for i < n + m, 



(3) for i < n + m, 



(4) for i < n + m — 1 , 



The estimates yield the following. 

Lemma 4.4. For any to of the form as in lemma \^~T\ and for any e > 0, there is an m > 
such that for any N and any m > m®, 



log 



Ctj{ft 
Cii(ft 



log 



^(ft 1 ) 



^ 2«— (n+m). 



log 



log 



^ 2*-( n + m )+i- 



a i (/3 JV+1 ) 



log 



< 



2«— (n+m)+l 
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Proof. The m the expression for converges, hence there is an ni\ > 1 such that 

the tail of the sum Yli>n+m a i a 2 ■ ■ ■ log ^ < f^. We will show how to choose m > mi 
to satisfy the conclusion of the lemma. 

We bound the influence of the change from (3 1 to f3 N using lemma 14.31 parts |2 and 0] 
The influence on each of the "head elements" (i < n + mi < n + m — 1) is bounded by 



log 



i-l 



^ ^ y ^—{n+m) _|_ 2*-(™+"i)+i <^ 2* _ ( n + m )+ 2 < 2 



mi +2— m 



^(/^...^(/^log^ 

By making m sufficiently large (i.e. by choosing a sufficiently large mo we can ensure that 

e a 1 (J3 N )...a i - 1 (J3 N )loz ! 



< 



16$(/31) ai (/F)... log 



1 



< 1 + 



hence 



16$(/?^ 
1 



< 



16$(^) / ' "a,^ 1 )' 

Adding the inequality for i = 1, 2, . . . , n + mi — 1 we obtain 

n+mi-l 1 n+mi— 1 -. 

x; ^...a^^iog^^- ^ ^(^...^(^log-^ 



8=1 



1=1 



16$(/3 



- £ a 1 (/5 1 )...a,_ 1 (/5 1 )log-i-< 



< 



n+mi-l 



i=l 



16' 



Hence the influence on the "head" of $ is bounded by 



To bound the influence on the "tail" we consider three kinds of terms 



a 1 (j3 N )...a i - 1 (j3 N )\o. 



(recall that i = n + m is not in $~). 
For n + mi < i < n + m — 2: 



: n + mi < i < n + m — 2, 2 = m + n — 1 and z > m + n + 1 



log 



ai^ 1 ) . . . ai-ii/3 1 ) log ^^ry 



a x (/?")... a*_ x (/3") log ^ 



i-l 



^ ^ ^ 2J~( TO + m ) _|_ 2*~( n + m )+ 1 < 2 i_ ( n+m ) +2 < 1 



Hence in this case each term can increase at most by a factor of e. 

For z = n + m — 1 Note that the change decreases log — - — so that log 

lQ g ^OT ' heilCeWeIiaVe 

a^.-.o^^bg^ a 1 (/^)...a i _ 1 03 JV ) 

log — 7^ TTrm 1 — ^ lo s — — < 



a„ +m _l(/3 Ar ) 



< 



ai^ 1 )... a,-!^ 1 ) log 
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n+m—2 

E2j~(n+m) < _ 
2' 

i=i 

Hence this term could increase by a factor of y/e at most. 

For i > n + m + 1: Note that Oj for j > n + m are not affected by the change, and the 
change decreases a n+m , so that a n+m (/3 N ) < a n+m {l3 1 ). Hence 

ai (n...«,- 1 (niog^ _ ai (^)...g w+m (^) 
° g aiOT • • • a^-it/? 1 ) log ^ ° g . . . a n+m (^) ~ 



So in this case each term could increase by a factor of e at most. 

We see that after the change each term of the tail could increase by a factor of e at most. 
The value of the tail remains positive in the interval (0, ||] , hence the change in the tail is 
bounded by f§ < f§- 

So the total change in $~ is bounded by 

change in the "head" + change in the "tail" < — + — = -. 
6 6 16 16 4 

□ 

Lemma 14.41 immediately yields: 
Lemma 4.5. For any e and for the same rjio(s) as in lemma \4~4[ for any m > itlq and N, 

|$-(^)-$-(^ +1 )|<|. 

Denote = ao(u)ai(u) . . . cn n+m _i(u/) log - — —rjr = $(u) — &~(u). We are now 

ready to prove the following. 

Lemma 4.6. For sufficiently large m, for any N, 

$ 1 (f3 N+1 )-$ 1 (f3 N )< £ -. 

Proof. According to lemma 14. HI part Q we have 

a l (/3»+ 1 )...a^- 1 {p»+ 1 ) " + "'- J 



log 



< ^ 2 *- {n+m) /iV < — . 



N 

8=1 



Hence . . . a n+rn ^(f3 N+1 ) < a^ N ) . . . a n+m _ 1 (/3 JV )e 1 / 7V , and 



lo S ^ + 0) 
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Hence 

l ° + iv' log(Af+l/0) J 

We make the following calcualtions. Denote x = ^^jv+t/^f ' t ' ien = A+l+1/0, 

(AT + Xj^y- 1 = 1 + ] r=i s < e^+W. N + 1/0 > e 1 / 3 , and so x - 1 < -j^rj-r < 4, thus 



JV+l/</> ^ ^ It 1 ^ ^ ^ iV+1/0 AT' 

3_ 

/V 

It is not hard to see that 0:^-10^ < 1/2 for all > 1, and we have 



1 /l 



(n+m-l)/2 



Thus 



$ 1 (/? 7V ) = a 1 (/3 JV )... a n+m _ a log -— < 1^) log(iV + l/0). 

- *-on < «■(/>") ( a + ^) lo g;;t / l // ) - 1 ) < 



-j \ (n+ro-l)/2 xj\(n+m-l)/2 -j. 



log(iV + 1/0) ((1 + e/AT)(l + 3/JV) - 1) < ( - ) log(iV + 1/0)-^ . 



Since ^ G o(l/ log(A/" + 1/0)), this expression can be always made less than | by choosing 



14 

A 7 " 

m large enough. □ 
Since $ = $~ + <3> 1 , summing the inequalities in Lemmas 14 . 51 and 14 . 61 yields the following. 

Lemma 4.7. For sufficiently large m, for any N, 

$(/5 7V+1 ) - $(/3 7V ) < e. 
It is immediate from the formula of <&(/3 N ) that: 

Lemma 4.8. 

lim $(/5 7V ) = 00. 

N— >oo 

We are now ready to prove lemma 14.11 

Proof, (of lemma l4~Tj) . Choose m large enough for lemma l4~7l to hold. Increase A" by one 
at a time starting with N — 1. We know that $(/3 1 ) = 3>(a?) < <&(u) + e, and by lemma 
Ol there exists an M with $(/3 M ) > $(co>) + e. Let A" be the smallest such M. Then 
$(^-1) < + £) anc i by lemma I4~7I 

$( / 9 JV ) < ${p N ~ x ) + e < + 2e. 

Hence 

+ £ < < $H + 2e. 

Choosing /3 = completes the proof. □ 
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The second part of the following Lemma follows by the same argument as Lemma 14.41 by 
taking N > 1 to be an arbitrary real number, not necessairily an integer. The first part is 
obvious, since the tail of u has only l's. 

Lemma 4.9. For an to = f3 l as above, for any e > there is an mo > ; such that for any 
m > mo, and for any tail I = [a n+m , a n+m+ i, . . .] if we denote 

(3 = [a±, ci2, . . . , a n , 1, 1, . . . , 1, a n+m , a n + m +i, . . .], 



then 



and 



ai(p 1 )a 2 (p 1 )...a i - l (p 1 )]og—^-r<e, 

t>n+m v ' 

n+m— 1 1 

aitf 1 ) ■ ■ ■ oa-iiP 1 ) log ~ MP 1 ) ■ ■ ■ log 



8=1 



< £. 



We can now prove lemma l4~2l 
Proof, (of lemma l4~2"|) . Applying lemma 1431 with | instead of e, we get 

$(/3 J ) - $(u;) = ^{"head"(/3 J ) - "head"(u;)} + "tail" (/3 J ) - "tail" (a;)} > 

-£{«taiF(w)} > — -- = -:. 



--V{"tairH}>---- 

□ 



We will need a computable version of Lemma 14.11 for modifying the conformal radius of 
the corresponding Julia set. 

Lemma 4.10. For any given initial segment I = (do, cii, . . . , a n ) and m-o > ; write u = 
[oq, at, . . . , a n , 1,1,1,...]. Then for any e > 0, we can uniformly compute m > itlq and an 
integer N such that if we write [3 = [a , Oi, . . . , a n , 1, 1, . . . , 1, N, 1, 1, . . .], where the N is 
located in the n + m-th position, we have 

(4.1) r(u) - 2e < r{/3) < r{u) - e, 
and 

(4.2) > 

Proof. We first show that such m and iV exist, and then give an algorithm to compute 
them. By Lemma 14.11 we can increase $(u;) by any controlled amound by modifying one 
term arbitrarily far in the expansion. 

By Theorem 12.101 / : 9 i— > + logr(0) extends to a continuous function. Hence for 
any Eq there is a 5 such that |/(x) — f(y)\ < £o whenever \x — y\ < 5. In particular, there 
is an mi such that \f(/3) — f(u)\ < £o whenever m>mi. 

This means that if we choose m large enough, a controlled increase of $ closely corre- 
sponds to a controlled drop of r by a corresponding amount, hence there are m > m and 
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iV such that (|4.1j) holds. (|4.2jl is satisfied almost automatically. The only problem is to 
computably find such m and N. 

To this end, we apply Lemma 13.11 It implies that for any specific m and N we can 
compute r(/3). This means that we can find the suitable m and N, by enumerating all 
the pairs (m, N) and exhaustively checking (|4.1|) and ()4.2|) for all of them. We know that 
eventually we will find a pair for which (j4.1j) and (|4.2|) hold. □ 

5. Proving the Main Theorem 

There are countably many oracle Turing Machines. Let us enumerate them in some 
arbitrary computable fashion Mf, M% , ... so that every machine appears infinitely many 
times in the enumeration. Recall that r{6) is the conformal radius of the Siegel disk 
associated with the polynomial Pg(z) = z 2 + e 2md z, or zero, if 9 is not a Brjuno number. 

We will argue by induction. On each iteration % of the argument we shall maintain an 
initial segment Ii = [ciq, a±, . . . , ajvj an interval Hi = [li, rj], and li = £{Hj) = — l{ such 
that the following properties are maintained: 

(5.1) Ti = r( 7i ), where 7; = [I h 1,1,.. .], 
and 

(5.2) For any (3 = [I h t Ni+1 , t Ni+2 , . . .] with r(j3) G [l h n], 

the machine Mf requires at least the time h(2 \— log£j]+l) to compute the y-approximation 
to J(P P ). And 

(5.3) for i > 1, $(/3) > $( 7i _!) - 2-( ? - 1 ), for any (3 = [I h t Ni+1 ,t Ni+2 , . . .]. 

Moreover, the intervals we construct are nested: [k,ri] C [k-i, Tj-i], and the sequence Ii 
contains as the initial segment. The numbers 2 \— log£{\ form a strictly increasing 
sequence. 

For the basis of induction, set Iq = [1], r = r(7 ) < 2 (by the Schwarz Lemma) and 
— r o/2, where 7 = [1, 1, 1, . . .]. Then for i = condition ()5.1|) holds by definition and 
conditions (|5.2|) and (|5.3|) hold because they are empty. 

The induction step. We now have the conditions (|5.1|) . (j5.2j) and ()5.3|) for some i and 
would like to extend them to i + 1. 

Consider the machine Mf +1 . Set = Simulate M^_ x on 7i for at most /i(2 [— log + 

1) steps to compute J 7i with precision -y 3 -. The machine reads at most /i(2 [— log£ i+ i] +1) 
bits of the input, and we can compute mo such that this run does not distinguish between 
7i and 7 = [I h 1, 1, . . . , 1, N mo +i, N mo+2 , . . .]. There are two cases: 

Case 1: does not terminate in the assigned time, or does not output a proper set. In 

this case, we proceed by setting Jj +1 = [Jj, !,...,!] (with l's up to position m ), 7i+i = 7j, 
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r i+ i = Tj, and U + \ = r i+ i — £i+±. By lemma we can choose sufficiently many l's in I i+ i, 
so that for any /3 beginning with Ij+i, we have > $(7i) — 2~\ 

Case 2: outputs a set S. Compute the conformal radius r(S). Considerations of 

Schwarz Lemma imply that for any quadratic Siegel disk, r(A) < 2. Using the above 
consideration to bound the constant in Lemma f2. 121 we know that for any Julia set J(P W ) 
which is £? +1 -accurately described by S, we have \r(u) — r(S)\ < Ay/2£ i+ i < 6£ i+ i. Again, 
there are two cases (if both hold, it doesn't matter which way to proceed): 
Subcase 2a: r^ — > r(S) + 8£ i+ i. In this case we proceed by setting J i+1 = [1$, 1, . . . , 1] 
(with l's up to position m ), 7 i+ i = 7i, r i+ i = n, and k +i = n+i - £i+i- 
Subcase 2b: U + 2£ i+ i < r(S) — 8£i + \. By Lemma 14.101 we can select = [Ii+i, 1, 1, . . .] 
by modifying ji at an arbitrarily far position, and set r i+ i = r(ji + i) so that $(7^+1) > $(7«), 
n+i >k + £i+i and [r i+ i - £ i+ x, r i+1 ] (7 [r(S) - 8£ i+ i, r(S) + 8£ i+ i] = 0. The number r m is 
computable since it is the conformal radius of a noble Siegel disk. Set k + \ = r»+i — 
We see that the induction is maintained for these parameters. 

In either subcase, by Lemma [4.21 we can add sufficiently many l's to Jj+i, so that for any 
(3 beginning with J i+1 , we have $(/3) > $(7i) — 2~ l , and condition (|5.3|) is satisfied. 

Lemma 5.1. Denote 7 = limi_ i . 00 7j. T/ien t/ie following equalities hold: 

$(7) = lim $(7i) and r(7) = lim r(7j). 

Proof. By the construction, the limit 7 = lim7j exists. We also know that the sequence 
r {li) = r i converges uniformly to some number r, and that the sequence $(7«) is monotone 
non-decreasing, and hence converges to a value ip (a priori we could have ip = 00). By the 
Caratheodory Kernel Theorem (see e.g. |Pomj ) . we have r^j) > r > 0, so ip < 00. On the 
other hand, by the property we have maintined through the construction, we know that 
$(7) > $(7i) - 2~* for all i. Hence $(7) > ip. 
From |BCj we know that 

(5.4) ijj + logr = lim($(7 i ) + logr(7j)) = $(7) + logr(7). 

Hence we must have ip = $(7), and r = r(j), which completes the proof. □ 

The conformal radius r(7) is computable, since the convergence r(ji) r(y) is uniform. 
Thus Jp i is also computable by Theorem 12.111 By construction, it satisfies all of the 
required properties. Note that the value 7 itself is also computable. 
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